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The XENONIOO and CRESST experiments will directly test the inelastic dark matter explanation 
for DAMA's 8.9cr anomaly. This article discusses how predictions for direct detection experiments 
depend on uncertainties in quenching factor measurements, the dark matter interaction with the 
Standard Model and the halo velocity distribution. When these uncertainties are accounted for, an 
order of magnitude variation is found in the number of expected events at CRESST and XENONIOO. 



Predictions for direct detection experiments require a 
wide-range of assumptions concerning the astrophysical 
properties of the dark matter, as well as its interactions 
with the Standard Model (SM). These theoretical un- 
certainties are compounded by additional experimental 
challenges that arise from the nature of low energy exper- 
iments. Ultimately, it is necessary to know the scattering 
rate for dark matter off SM nuclei in detectors. There 
are many unknown physical quantities that go into this 
prediction and they are often benchmarked to values in 
specific studies. However, in light of a potential signal, 
the verification process requires a more systematic study 
of these unknowns in order to have a complete picture of 
the range of consistent theories. 

Inelastic Dark Matter (iDM) serves as a case study 
for this new treatment of uncertainties and shows how 
marginalizing over astro and particle physics quantities 
leads to at least an order of magnitude variation in detec- 
tion prospects at upcoming experiments. Inelastic dark 
matter is an elegant explanation for DAMA's on-going 
8.9(7 annual modulation signal [T], resolving the inconsis- 
tency of this signal with the plethora of null direct de- 
tection experiments [2H3 ■ In the iDM framework, a dark 
sector particle up-scatters off the detector's target nu- 
cleus to a higher mass state [S]. To explain the DAMA 
anomaly, an 0(100 keV) mass splitting is required for 
weak-scale dark matter. 

IDM requires a minimum velocity to up-scatter to the 
more massive state, which depends on the mass of the 
target nucleus, m^, the reduced mass of the nucleus-dark 
matter system, /xat, the mass splitting, 5, and the recoil 
energy, Er^ of the nucleus: 
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VE{t) accounts for the boost to Earth's rest frame |10| . 
The differential scattering rate is larger for heavier target 
nuclei because Vniin is reduced. 

The spin-independent differential cross section can be 
parameterized as 
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where a^ is the dark matter-nucleon cross section at zero 
momentum transfer, ^n is the dark matter-nucleon re- 
duced mass, and q^ = 2mNEfi is the momentum transfer. 
The constants /p^n parameterize the coupling to the pro- 
ton and neutron, respectively, and are set to /p — fn — ^ 
throughout. The dependence of the cross section on the 
nuclear recoil energy comes from the dark matter and 
nuclear form factors, Fdml?^) and F^iq^)- F^i^{q^) de- 
scribes non-trivial behavior at low momentum transfer in 
models where higher dimensional operators contribute to 
the scattering [TTHI3]- F^{q'^) is the Helm/Lewin-Smith 
nuclear form factor [9]. Analytic approximations to the 
nuclear form factor can have substantial errors, particu- 
larly for heavy nuclei such as ^*'*W. The Helm/Lewin- 
Smith form factor is better behaved than other Helm 
parameterizations, but can still give errors of 25% for 
is^'W in the range Er = 10-40 keVnr [H]. Around 100 
keVnr, these errors can be as large as 60%. The impact 
of nuclear form factor uncertainties on predictions for di- 
rect detection has been addressed in the literature |15j . 
This work explores other sources of uncertainty, and for 
the rest of this paper the Helm/Lewin-Smith form factor 
is adopted. 

Typically, the vdf is taken to be Gaussian, isothermal 
and isotropic in the galactic frame. This 'Standard Halo 
Model' (SHM) is parameterized as 
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where f{v) is the local velocity distribution function 
(vdf) for the dark matter halo in the galactic frame, and 



where Wosc is the galactic escape velocity and vq is the 
velocity dispersion. The range of escape velocities is con- 
strained by the RAVE stellar survey: 480 < Wesc £ 650 
km/s |16j . and no constraints are placed on vq. The stan- 
dard procedure when evaluating direct detection rates is 



to assume a SHM distribution with benchmarked values 
for vq and Uesc- The soUd blue curve in Fig. [l] shows the 
expected tungsten recoil spectrum for this vdf; the bulk 
of events occur between 10-40 keVnr. 

While previous studies have looked at the effect of 
varying vq and Wosc within the SHM [T71 [TS] , none have 
fully marginalized over both dark matter and halo profile 
uncertainties. In addition, numerical N-body simulations 
indicate significant departure of the vdf from the SHM 
hypothesis [TM2^ . especially in the high velocity tail. 
Because very little is known about either the vdf or the 
dark matter model, experimental analyses should be de- 
signed to cover a wide range of possibilities. In this paper, 
a scan over the parameter space for iDM is performed, 
where we marginalize over the dark matter parameters 
{m,d,a) and halo velocity parameters (I'oj^'esc^'strcam), 
and set constraints by a global x^ analysis ^5].^ 

The predicted number of events at CDMS [2], 
ZEPLIN-H jg, ZEPLIN-III [4], CRESST-H [5], 
XENONIO 0, EDELWEISS [7], and the XENONIOO 
calibration run |26| are included in the x^ as well as the 
annual modulation amplitude in the first twelve bins of 
DAMA (2-8 keVee) [I]. The high energy bins from 8-14 
keVee are combined into a single bin with modulation am- 
plitude -0.0002±0.0014 cpd/kg/kcVce. Any model that 
over-predicts the number of events at the null experi- 
ments by 2(7 is excluded. 



Inelastic Dark Matter at CRESST 

The standard assumption is that the iDM interpreta- 
tion of DAMA's signal can be confirmed or refuted by 
any experiment with a target nucleus heavy enough to 
cause the inelastic transition [25]. For DAMA, the in- 
elastic transition occurs through scattering off the ^^^I 
nucleus in the Nal(Tl) target. Naively, any experiment 
with a target mass greater than ^^^I could provide a suf- 
ficient test. Two upcoming experiments fall into this cat- 
egory: XENONIOO i26j and CRESST ^^, which use 
"iXe and ^^^^W, respectively. XENONIOO has currently 
released results from a calibration run of 11.2 live days 
during Oct-Nov 2009. They report 161 kg-d effective 
exposure and have observed no events in their accep- 
tance region between 4.5 - 40 keVnr. CRESST, which 
consists of nine detectors of CaW04 and one detector 
of ZnW04, has shown preliminary results in the energy 
window from 10 - 40 keVnr obtained from summer 2009 
until the present; however, the exposure was not reported 
[28]. 
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^ In some circumstances, maximum gap techniques provide tighter 
limits than Poisson statistics for null experiments |15II24) . How- 
ever, Poisson statistics are used in this paper due to the complex- 
ity of combining a x^ for DAMA with multiple max-gap tests to 
get a global limit. 



FIG. 1: CRESST spectra for: regular iDM (blue), FFiDM 
with Fdm oc En (green), DM stream (red) for Qi — 0.085. 
Inset: Qi = 0.07 (dashed), Qi = 0.06 (dotted). 



CRESST provides a unique experimental environment 
for testing iDM because it has the heaviest target nu- 
cleus of all current direct detection experiments. ^®^W is 
expected to be highly sensitive to inelastic scattering be- 
cause its velocity threshold is a factor -y/mi/mw ~ 0.83 
lower than iodine. As a result, a larger fraction of the 
halo can up-scatter off of ^®^W and one would expect a 
larger scattering rate compared to lighter targets. How- 
ever, an additional complication arises due to the large 
radius i?w of a ^*'*W nucleus. In particular, when the mo- 
mentum transfer is q ^ l/i?Wj the dark matter probes 
the size of the nucleus and the scattering is no longer 
coherent. Therefore, the scattering rate is suppressed at 
recoil energies Er oc l/(2mNi?w)- This suppression oc- 
curs at lower recoil energies for -'^^''W, as compared to 
^^^Xe (55 keVnr versus 90 keVnr, respectively). 

The fact that the first zero of the ^*^W form factor 
occurs at such a low recoil energy highlights an impor- 
tant challenge for CRESST. The typical recoil energy an 
iDM particle deposits on a ^^''W nucleus is 0{5 fj,/m^^) ~ 
0(50 keVnr), where the signal is suppressed due to loss of 
coherence. The form factor suppression is evident in the 
recoil spectra of Fig. [I] which also illustrates the effects 
of three additional sources of uncertainties: DAMA's en- 
ergy calibration, the dark matter interaction with the 
SM, and the velocity distribution profile. Variations in 
any of these three sources can significantly alter iDM's 
rate at CRESST and ultimately affect the final exposure 
that will be required for CRESST to exclude the DAMA 
iDM hypothesis. 



DAMA 's Energy Calibration 

DAMA only detects its nuclear recoil events with scin- 
tillation light. Nuclear recoils typically deposit only a 
small fraction of their energy into scintillation. The 
quenching factor for ^^^I , Qi, relates the measured elec- 
tron equivalent energy (given in keVee) to the nuclear 



recoil energy (given in keVnr): 
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where the energy dependence of the quenching factor is 
left explicit. Most studies assume a constant quench- 
ing factor for iodine from '^ 10 — 100 keVnr, with the 
standard value taken to be Qi — 0.085. However, 
there are large experimental uncertainties in measure- 
ments of Qi [29 . The four primary ones [30H33] give 
0.05 < Qi < O.io. The study in [32] gives the small- 
est error, however its measurements are calibrated with 
60 keV gamma rays, in contrast to the 3.2 keV electrons 
that DAMA uses [30; . This difference reduces the central 
value of [32^ by roughly 10% and induces larger system- 
atic effects. 

Lowering iodine's quenching factor effectively shifts 
DAMA's signal to higher nuclear recoil energies, favor- 
ing slightly larger values for the iDM mass splitting 
(100 < 5 < 180 keVnr for Qi = 0.06). Consequently, 
the predicted signal at other experiments is also shifted 
to higher nuclear recoil energies. In addition, the spec- 
tral shape is broadened, because DAMA's reported rate 
is in units of cpd/kg/keVee. 

Fig.[l]illustrates how the ^*^W recoil spectrum changes 
as Qi is reduced from 0.085 to 0.06, and Fig. [2]shows the 
average annual rate for CRESST's low and high energy 
range, assuming a SHM profile. The shift of the iDM sig- 
nal to higher recoils translates into a significant reduction 
in CRESST's average annual rate in the low recoil win- 
dow of 10 — 40 keVnr and a substantial enhancement in 
its rate in the high recoil range from 40 — 100 keVnr. 

The effect of marginalizing over both the particle and 
halo parameters is substantial. As an illustration, if only 
the mass and cross section are marginalized over (with 
fixed S = 120 keV, wq = 220 km/s, Wosc = 550 km/s), then 
the average counts at CRESST per 100 kg-d range from 
~ 17 — 36 in the lower energy window and ~ 0.3 — 0.8 in 
the high energy window for Qj = 0.085. At XENONIOO, 
the predicted number of events (per 1000 kg-d) ranges 
from '-^ 33 — 100, covering only a subset of all allowed 
possibilities illustrated in Fig. [3j 

Dark Matter Interaction 

The identity of iDM is unknown and its interactions 
with the SM may not occur through renormalizable op- 
erators. Non-renormalizable operators typically result in 
matrix elements with non-trivial dependence on the mo- 
mentum transfer q. These can be parameterized by an 
effective DM form factor [TMll [H [34H36] , 
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where g^ = E^, \q\ = -/2mjv^B/{, and A is an arbi- 
trary mass scale. Standard iDM assumes that the con- 
stant n, m = 0, term dominates the expansion. Models 
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FIG. 2: Average counts at CRESST per 100 kg-d for regular 
iDM (blue), FFiDM with F^miq) oc Er (green), and DM 
streams (red) . The effect of lowering the quenching factor is 
illustrated for Qi — 0.07 (dashed blue) and Qi = 0.06 (dotted 
blue). The contours enclose all points with x^ ^ 18- 



that have an interaction mediator with mass lighter than 
Odffl) are dominated by co,-2- Composite iDM models 
have co,i ^ [571139) . Form factors that are dominantly 
n ^ can be realized through dipole or other tensor 
interactions [551 [3S] . 

Standard iDM (i.e., n,m = 0,0) and models with 
n =: 0, m 7^ have comparable rates at CRESST because 
the ratio of predicted events between these two scenarios 
scales as iVo,m/^o,o - (mwEw pcak/miEi poak)^™ ^ 1- 
DAMA's spectrum peaks at Ei peak — 35 keV while the 
tungsten spectrum at CRESST peaks at i?w peak — 
25keV. In contrast, interactions with n ^ 0,m = pre- 
dict substantially smaller rates at CRESST: A'^„^o/A'^o,o — 
(-E'w ponk/Ei poak)^" — (0.5)". This effect is illustrated in 
Fig.[|for n,TO = 1,0. 



Dark Matter Velocity Distribution 

There is little direct observational evidence for the 
DM density profile, and the velocity distribution is 
highly uncertain. While most studies assume a Maxwell- 
Boltzmann vdf Q , N-body simulations indicate that this 
ansatz does not adequately parameterize the vdf |22j . 
The iDM spectrum is particularly sensitive to changes 
in the tail of the velocity distribution profile, which can 
arise from velocity anisotropics or from dark matter sub- 
structure that has recently fallen into the galaxy fT9j |40] . 
A vdf in which the high velocity tail is dominated by a 
stream of dark matter illustrates how changes in the local 
vdf alter iDM predictions. This scenario can significantly 
lower the number of expected events; other possibilities 
for the velocity profile will result in numbers of events 
between the SHM and stream expectations. 

Streams of dark matter are characterized by low veloc- 
ity dispersion [19j . Here, streams will be parameterized 
as dispersionless vdfs that have an arbitrary incident an- 
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FIG. 3: Average counts at CRESST (per 100 kg-d) versus 
XENONIOO (per 1000 kg-d) for regular iDM (blue), FFiDM 
with Fdmiq) oc Er (green), and DM streams (red). The effect 
of lowering the quenching factor is illustrated for Qi = 0.07 
(dashed blue) and Qi = 0.06 (dotted blue). The contours 
enclose all points with x^ ^ 18- 



gle. The distribution profile is f{v) — 5^{v — ^stream), 
with V and ^stream given in the frame of the sun. The dif- 
ferential scattering rate is obtained after boosting to the 
Earth's frame, and depends on the recoil energy through 
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where VE{t) is the Earth's velocity in the frame of the so- 
lar system. The rate would be constant if not for the nu- 
clear form factor that shapes the distribution and yields 
a highly peaked spectrum as illustrated for tungsten in 
Fig.[l] 

The stream's velocity and incident angle relative to the 
Earth are marginalized over and constraints on the phase 
and higher harmonics of the annual modulated rate are 
applied. These are set by the spectral decomposition 
of DAMA's modulated rate, which restricts the modu- 
lation to peak at May 24**^ ± 7.5 days and constrains 
the power spectrum at a frequency w = 2 yr~^ to be 
P(2 yr~^) < 0.05P(yr~i) [1,. We apply these additional 
cuts on the models that survive the x^ minimization. 

Fig. [2] shows the iDM predictions for CRESST when 
the tail of the vdf is dominated by a dispersionless 
DM stream. When the dark matter stream is nearly 
head-on in the summer, very few events are expected 
in the fall and winter months, which is consistent with 
the late-year running of the XENONIOO calibration run. 
Marginalizing over stream parameters shows that veloc- 
ities Wstroam ~ 400 km/s are favored, and that the DM 
incident angle with respect to the velocity of the Earth 
on June 2 can be as large as 75°, although 90% of the 
consistent models have incident angle ^in < 50° and 75% 
have 9in < 36°. Fig.|2]shows that the annual average rate 
at CRESST for dark matter streams can deviate dramat- 
ically from the SHM case and highlights the importance 
of marginalizing over parameters and considering differ- 



ent velocity profiles when making predictions for direct 
detection experiments. 



XENONIOO Prospects 

The previous section discussed how uncertainties in 
energy calibration, dark matter interactions, and vdfs 
significantly affect the range of predicted events at 
CRESST. For 0(100 kg-d) exposures, the number of 
events in the low energy window can vary over an or- 
der of magnitude from 3 — 30, while the number in the 
high energy window ranges from 0.1 — 10. It is evident 
that having significant exposure over the full nuclear re- 
coil band where the iDM signal is expected to dominate 
(10-100 keVnr) is essential for refuting or confirming a 
potential signal. 

The XENONIOO experiment will accumulate large ex- 
posures 0(3000 kg-d) in their current data run. Com- 
pared to CRESST, it has better coverage of the relevant 
nuclear recoil band because the ^"^^Xe target is not af- 
fected by form factor suppression at energies below 90 
keVnr. The dominant variation in predictions for spin- 
independent iDM scattering rates in ^'^^Xe arises because 
DAMA's modulation fraction is not known and only con- 
strained to be greater than 2% [T]. Larger modulation 
fractions at DAMA imply a proportionately smaller rate 
at XENONIOO. 

In their calibration run, XENONIOO demonstrated the 
potential to be a "zero background" experiment. The 
fact that no events were seen in their acceptance region 
implies a lower bound on the modulation fraction for iDM 
of 0(40%). Their next data release, which will include 
summer data, will directly test the iDM parameter space 
still allowed. Fig. [3] illustrates the average number of 
expected events for XENONIOO (per 1000 kg-d), includ- 
ing the uncertainties in Qi, the DM form factor, and the 
vdf. Again, there is an order of magnitude uncertainty, 
with the predictions ranging from 20-200 counts per 1000 
kg-d. For the expected exposure of XENONlOO's data re- 
lease, the minimum number of events predicted by iDM is 
~ 60, which will be enough to confirm or refute the spin- 
independent iDM scenario in a conclusive and model- 
independent manner. 



Conclusions 

The process of testing the DAMA anomaly highlights 
many of the challenges inherent to direct detection ex- 
periments. In addition to determining the properties of 
the unknown dark matter particle, direct detection ex- 
periments must also consider the unknown flux of the 
incident dark matter, as well as uncertainties in convert- 
ing a signal from one target nucleus to another. 

The predictions for both the CRESST 2009 run and 
XENONIOO 2010 run show an order of magnitude uncer- 
tainty. The nuclear form factor for ^^^W, when combined 



with additional theoretical and experimental uncertain- 
ties, will likely prevent CRESST from refuting the iDM 
hypothesis with an exposure of 0(100 kg-d) in a model- 
independent manner. XENONIOO, on the other hand, 
will be able to make a definitive statement about a spin- 
independent, inelastically scattering dark matter candi- 
date. StiU, the CRESST 2009 data can potentially con- 
firm iDM for a large range of parameter space. In case of 
a positive signal, the combined data from CRESST and 
XENONIOO will start probing the properties of the Milky 
Way DM profile and the interaction of the SM with the 
dark matter. 
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